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Abstract 

The numerical solutions of the non-linear evolution equation are shown to display the 
"geometric" scaling recently discovered in the experimental data. The phenomena hold both 
for proton and nucleus targets for all x below 10~^ and 0.25 GeV^ < < 2.5 x 10^ GeV^. 
The scaling is practically exact (few percent error) in the saturation region. In addition, an 
approximate scaling is found in the validity domain of the linear evolution where it holds 
with about 10% accuracy. 

Basing on the scaling phenomena we determine the saturation scale Qs{x) and study 
both its x-dependence and the atomic number dependence for the nuclei. 
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1 Introduction 



The experimental data on the structure function F2 were recently discovered to display exciting 
phenomena called the "geometric" scaling Namely, the total 7*p cross section is not a function 
of two independent variables x and Q but rather a function of a single variable r = Q/Qs{x). 
The function Qs{x) is a new scale called the saturation scale. The scaling holds experimentally 
with 10% accuracy in the whole kinematic region of x < 10~^ |jl|. 

These remarkable phenomena require theoretical explanation. In fact, the scaling behavior is ac- 
tually anticipated in high density QCD and is strongly related to the appearance of the saturation 
scale 0, ^ H, |]. During the interaction a parton cascade is developed. When the parton density 
(we mean the packing factor) is not large the transverse momenta of the partons are strongly 
ordered. Such a system evolves according to the linear DGLAP equation which describes the 
gluon emission. As a result of this radiation the number of partons rapidly increases. However 
in the high parton density phase annihilation processes become significant and they suppress the 
gluon radiation resulting in the saturation of the density. This scenario happens at the saturation 
scale Qs{x) [0, 0, which has a meaning of an average transverse momentum of partons in the 
cascade. At photon virtualities below Qs{x) the ordering in the momenta does not persist any 
more and all partons have the same momenta Qs{x). This is a domain where the evolution cannot 
be described by a linear equation. A non-linear evolution should be used instead. The scaling is 
a property of this kinematic region. It just says that the very same average momentum Qs{x) in 
the cascade can be approached from two directions either by varying the virtuality Q at fixed x 
or vise versa. The scaling phenomena are expected aX Q < Qs{x). Furthermore, the saturation 
scale Qs{x) can be defined as a scale at which the scaling breaks down. 

As a result of the above discussion we conclude that the scaling phenomena should be a con- 
sequence of the non-linear evolution whith the non-linear effects switching on at the saturation 
scale. Numerous theoretical efforts to understand theoretically the mechanisms responsible for 
the parton saturation |P, |^, |10|, |Tl|, |12|, |13| lead finally to the very same nonlinear evolution 
equation |]^, ||, |lT], 0, |1^, [1^, [15[ . This equation has a most transparent form in the color dipole 
approach: 



iV(xoi,r;6) 



iV(xoi,lo;&) exp 
Cpas 



Y 



dyexp 




+ 



X 



(1.1) 



^01 



^02 ^12 



2 A^(xo2, b - ^xi2) - A^(xo2, b - ^xi2)iV(xi2, y; b - ^xo2; 



r^,x; b) whith a'^^poie being the elastic ampli- 



The equation is written for N{r±, x; b) 
tude for a dipole of size scattered at the impact parameter b. The rapidity Y = —Inx and 
Yo = — InxQ. The ultraviolet cutoff p is needed to regularize the integral, but it does not appear 
in physical quantities. In the large Nc limit (number of colors) Cp = Nc/'i. 
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The equation (|1 . 1| ) describes the following physical picture. The evolution kernel 
probability for the dipole of size xio to decay into two dipoles of sizes Xi2 and xo2- Then these 
two dipoles interact independently with the target (linear term in the equation). The nonlinear 
term in the evolution takes into account the Glauber corrections for the interaction. These 
corrections are due to the screening between the two dipoles, and hence they contribute with the 
negative sign. 

It can be seen from the form of the equation (|1 . 1| ) that it contains no information about target. 
The only dependence on a target is coded in the initial conditions of the evolution at some initial 
value xq. We take for the initial conditions the Glauber-Mueller (GM) formula, which is proven 
to be correct initial distribution for nuclear targets [1^. For the proton target we also use the 
Glauber - Mueller formula. However, in this case the procedure is less justified theoretically and 
we rely in our choice on the fact that this formula describes well the experimental data at not 
too low X. The initial conditions are: 



A^(xoi,xo;&) = A^GM(xoi,a;o;6) 



with 



iV'GAf(xoi,a;; 6) = 1 - exp 



(x,4/x^i)5(b) 



;i.2) 



;i.3) 



The equation ( |1.3| ) takes into account the multiple dipole-target interaction in the eikonal ap- 
proximation |jl6|, |T^, |T^. The gluon density xG^'^^^^ is a solution of the DGLAP equation [||. 
The function S{b) is a dipole profile function inside the target, while R stands for its radius. 

Solutions to the equation ( p,.l| ) were studied in asymptotic limits in Refs. |]. Numerical 
solutions of the equation (|1.1|) were reported in the Refs. ||15|, |19|, We continue studying 

the properties of the solutions obtained in the Refs. [|1^, |21[ . In the present work we concentrate 
on the scaling phenomena displayed by the solutions of (|1 . 1| ) . Indeed, it was shown in the Ref. 



that in the double logarithmic approximation, the solutions of the equation ( |1 . 1| ) scale with 



a good accuracy in a wide high energy region. The recent paper pO[ reports on the numerical 
observation of the scaling phenomena. 

The main goal of the present paper is to discover the scaling phenomena in the solutions for 
both proton and nuclei obtained in the Refs. |21[. Basing on this property we determine the 
saturation scale Qs{x) and for nuclei we study its A dependence. 

The paper is organized as follows. The next section (2) is devoted to the scaling phenomena 
and saturation scale for the proton. The scaling on nuclei is investigated in section 3. The final 
section (4) concludes the work. 



2 Scaling and saturation scale for proton 



In the recent paper |T9[ the nonlinear evolution equation ( |1.1| ) was solved numerically for constant 



value of the strong coupling constant as = 0.25. The goals of the present research are in further 
study of the physical properties displayed by the solutions of ( |1.1| ). 
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Compared to the Ref. [|19| in the present paper we shghtly modified the large distance behavior 
of the initial conditions at x = 10~^. In fact, no information about large distances is known. 
The GRV parameterization p2| entering the initial conditions ends up at distances ^ 0.5 fm. 



In the previous paper |]I9| the extrapolation to larger distances was done by a constant, which 
does not approach unity at the very large distances. Moreover, such initial conditions cannot 
be consistent with a scaling being a purely dynamical property of the evolution equation. There 
exists a transition region below x = Xq = 10~^ where the solutions of (|1.1| ) are sensitive to the 
initial conditions. In this transition region we do not expect to observe any scaling phenomena. 
The transition region is estimated to end up at a; ~ 10~^. Below x ~ 10"'^ the initial conditions 
are forgotten and the dynamics is governed by the pure evolution. Then the scaling sets in and 
indeed it is seen in the solutions obtained in the Ref. 



In order to eliminate the transition region, in the present work we extrapolate the large distance 
behavior of the initial conditions consistently with the asymptotics and the scaling. The following 
procedure can be suggested. We take a solution of Ref. |T^ at some x well bellow the transition 
region, say at x = 10~^. Appropriately rescaled this solution is used for the large distance 
extrapolation of the initial conditions at a; = 10^^. The above described improvement of the 
initial conditions modifies slightly the large distance behavior of the solution in the transition 
region 10~^ < a; < 10~^ and restores the scaling in this region. Note that in a sense our procedure 
already implies the scaling at very large distances. 

The solutions obtained N{r±,x) = N{r±,x]b = 0) are displayed on the Fig. (|l]). The different 
curves correspond to the different values oi x. As can be seen from the Fig. at any fixed x 




Figure 1: Solutions of the equation as a function 

of distance. The different curves correspond 
to solutions at X = 10^"^ (the upper curve), 
10-6 and so on down to x = 10 ^ (the lowest 
curve ) 
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the solution N behaves in a step like manner as a function of distance: at small distances it tends 
to zero, while at larger distances the saturation value unity is approached. 



2.1 Scaling phenomena 

In this section we study a possible scaling behavior of the solution A^. As has been mentioned, 
the double logarithmic approximation of the solutions of the master equation ( |1.1|) ||^ as well as 
general analyses of similar non-linear equations [0] predict this new scaling phenomenon in the 
saturation region r± > 1/Qs{x). In the saturation region this scaling implies the dipole-target 
amplitude to be a function of only one variable r = ■ Ql{x): 

N{r±,x) = N{t). (2.4) 
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Indeed, with proper rescaling of the variable r±_ all the curves in the Fig. (|T]) can be mapped one 
onto another. This is a manifestation of the scaling property ( |2.4| ). 

A rigorous numerical procedure for the scaling detection can be defined. It is useful to introduce 
the rapidity variable y = Inl/x. Let us define the following derivative functions assuming the 
scaling behavior ( p.4|) : 

Nrir^,x)^rl^^ = ^rlQ^x). (2.6) 

From the equations (|2.5|) and ( |2.6| ) one can see that the ratio Ny/Nr is a function of only one 
variable x: 

Ny dlnQ^Jx) 
Ra{r^,x) = ^ = — TT^- 2.7 

Our goal is to investigate the above ratio from the solutions obtained. In the Fig. (|^) the derivative 
functions Nr and Ny are plotted versus for various values of x. One can clearly observe 
similarity in the behavior of these functions. This is actually a sign of the scaling phenomena. 
Both functions iVr and Ny possess extremum points at which the derivatives with respect to r_L 



vanish. If the scaling behavior takes place then it follows from (|2.4| ) that both Nr and Ny are 
extreme at the same points. In fact, this condition is clearly observed with a very good accuracy 
(Fig. I). 
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Figure 2: The derivative functions Nr (dashed line) and Ny (solid line) as functions of the distance 
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In order to establish the scahng phenomena numerically we have to check if the function Ra is 
indeed independent. However, it is clear that we cannot expect exact numerical independence. 
So, a numerical criteria for the scaling existence has to be defined. In this paper we study scaling 
within the distance interval 0.04 GeV"^ <r±<10 GeV~^ that corresponds to 0.25 GeV^ < < 
2.5 X lO^GeV^. Since the experimental accuracy for the scaling is about 10% we define the 
following condition for its acceptance: 



max ARa 



max AA^, 



< 10%, 



(2.8) 



r,y 



where max ARa is a maximal variation (in percents) within the interval of interest of the function 
Ra with the distance r± at fixed x. The functions max ANr^y are similarly defined. The condition 
( |2.8| ) says that we accept for the scaling some small r± dependence of the function Ra in a scale 
of large variations of the functions Nr and Ny. 

The figure (^ presents the main results. The three lines correspond to functions Nr and Ny 
divided by their minimal values within the interval, and the function Ra multiplied by the factor 
20 to be seen on the scale. 
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Figure 3: The scaling as a function of the distance r±. The positive curves are Nr /Nvmin (dashed line) 
and Ny /Numin (solid line). The dotted line is 20 x Ra. 



The function Ra is clearly observed to be a very slowly varying function of r± for all values 
of X and r^. Though at fixed x the function Ra cannot be claimed to be exact constant, its 
variations with are very much suppressed comparing to the variations of the functions Nr and 
Ny. For example, at x = 10~^ within the given interval the function Ra changes by maximum 
15%. Within the very same interval both functions Nr and Ny change in several times. Then the 
relative fluctuation is much less than 10%, which according to the condition ( |2.8| ) confirms the 
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scaling. The phenomenon holds with a few percent accuracy and it improves at smaller x ~ 10~^ 
and in the deep saturation region. However to observe this scaling behavior in these regions is 
numerically more problematic since both derivatives Nr and Ny tend to zero. 

As was discussed in the Introduction, we should expect the scaling violation at distances of order 
1/Qs{x). The shorter distances are in the realm of applicability of the linear equation which is 
not supposed to display any scaling phenomena. Nevertheless, we observe the scaling actually to 
exist also at distances which are much shorter than the saturation scale. The above statements 
seem to contradict each other. We believe, however, that the resolution of the paradox is in 
the linear equation which in fact exhibits the approximate scaling behavior Unfortunately, 
this numerical coincident prevents us from determination of the saturation scale as a scale of the 
scaling violation. 



2.2 Saturation scale 

No exact mathematical definition of the saturation scale is known so far. In the Ref. ||T9|1 two 



definitions of the saturation scale were proposed and the solutions obtained from the equation ( |1.1|) 
were used for its determination. For the step-like function it is natural to define the saturation 
scale as a position where = 1/2: 

• Definition 1: 

N{2/Qs,x) = 1/2. (2.9) 
An alternative definition of the saturation scale is: 

• Definition 2: 

K = - ln[l - N{2/Qs,x)] = 1/2. (2.10) 

The latter definition is related to the 6-dependence of the solution and is motivated by the GM 
formula with n being the gluon packing factor[|. This definition is equivalent to N{2/Qs, x) ~ 0.4 
which predicts a somewhat larger saturation scale Qs{x) comparing with ( |2.9|) . The saturation 
scales obtained through the equations (|2.9| ) and ( ^.10[ ) are plotted in the figure (^a). 

The saturation scale can be deduced directly from the scaling property ( p.4|) which have been 
established. To this goal we can regard the equation (p.7|) as a new definition of the saturation 
scale: 

• Definition 3: 

— ' = Ra{r^,x) = const{r^) (2.11) 



^In the present paper as well as in the previous papers |l9[ ^ we do not deduce the saturation scale Qs 
but rather the dipole saturation radius Rg. The equality Qs = 2/Rs is motivated by the double logarithmic 



approximation. Though for the equation (1.1) this approximation is not justified, we still believe it to make 



reliable estimates provided Qs is sufficiently large. 
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Figure 4: The saturation scale Qg is plotted as a function of x. (a) - the scales obtained from the 
equations ( \2.{\ ) (solid line) and ( 2.1(\ ) (dashed line), (b) - the equation ( 2.1^ ) is used to 
determine the scale, (c) - the result obtained from the equation {2. IS). 



This definition allows us to determine the energy dependence of the saturation scale. Note from 
the Fig. (^) that the function Ra is practically independent of x, Ra ^ —0.7. Hence we obtain 
from (|2.11|) the power law dependence of the saturation scale on x: 

Qs{x) = Qsox-'i = g.oe-^'""; q = 0.35 ± 0.04. (2.12) 

The obtained dependence of the saturation scale on x is somewhat weaker than both the double 
logarithmic prediction q = 2asNc/ tt and the numerical result of the Ref. |]2D| but significantly 
stronger than the GW saturation model qcw = 0.288/2 ||23|. Unfortunately, the parameter Qso 
cannot be deduced from the scaling analysis only. In order to make some estimates we choose 
two reasonable values for Qso just fixing the saturation scale at x = 10~^: (5^(10"^) = 1 GeV and 
Qs(10~^) = 1.5 GeV. The obtained results are plotted in the Fig. (||,b). 

The physical meaning of the scale at which both Nr and Ny are extreme is that at this scale 
(which is a function of x) the nonlinear effects responsible for the saturation set in. Thus it is 
natural to suggest yet another definition of the saturation scale 



Definition 4: 



(9r^ dx 



ri=4/Q2(x) 




(2.13) 



It is important to stress that the definitions ( p.llD and (|2.13| ) are not equivalent, and no one is 
a consequence of another. However, if we suppose that k ~ ('"D^ in analogy with the GM 



formula then both definitions (|2.11|) and (|2.13|) are equivalent 



The Fig. (^,c) shows the saturation scale Qs{x) obtained from ( |2.13| ) as a function of x. The 
values presented are deduced with few percent errors. The Fig. (|^,c) predicts smaller saturation 
scales compared with the ones from the Fig. (^a). This fact can be naturally explained if we 
suppose again that n ~ {f±Y~'' ■ Then it is easy to show that the definition ( p.l3| ) corresponds to 
the condition k{x,2/Qs) = 1, which implies N{2/Qs,x) = 1 — 1/e ^ 0.63. Hence, the saturation 
is obtained at larger distances. 
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3 Scaling phenomena for nuclei 



Solutions of the equation ( |1.1| ) were obtained for nucleus targets in the recent paper All 
details about the solutions for six nuclei Auiqt, Ndi^Q, Moioo, ZnjQ, Ca^o, and Ne2o can be 
found there. In the present work the only modification we perform is again concerns the large 
distance extrapolation of the initial conditions. To this goal we use the Glauber formula for the 
initial conditions. For the nucleon cross section we use the result of the previous section. This 
way we discover that the Glauber formula extrapolates the large distance behavior of the initial 
conditions consistently with both the asymptotics and scaling. The solutions obtained NAir±,x) 
display similar step-like behavior as in the proton case (Fig. |1|). 

In this section we investigate the scaling phenomena for nuclei following the very same strategy 
as presented above. We start from computations of the derivative functions NrA and NyA for 
the nuclei, where the subscript A stands for the atomic number of a nucleus. In fact, the depen- 
dences quite similar to the Figs. (H) and (||) are obtained. The Fig. (H) presents an example of 
calculations for the most heavy nucleus Au. Since in our approach the solution for Ne is almost 



identical to the proton (see the Ref. |^ for the discussion) Ne nucleus displays exactly the same 
scaling phenomena as proton. 
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Figure 5: The scaling as a function of distance the r±). The positive curves are NrAu/NrAm 
NyAu/NyAumin- The dotted line is 20 x Rqau- 



and 



The scaling on gold is observed form the Fig. (^. The other nuclei display the very same 
phenomenon. Though for nuclei the numerical fluctuations are larger it holds within at least 10% 
accuracy with respect to the condition (|2.8| ). 
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Nuclei \ X 


10-^ 


10-^ 


10-5 


10-^ 


10-3 


Light 


0.20 


0.23 


0.26 


0.26 


0.29 


Heavy 


0.18 


0.18 


0.22 


0.22 


0.26 


All 


0.19 


0.20 


0.23 


0.23 


0.27 



Table 1: The power P2{x) for various values of 

X. 



The ratio RaA{x) is almost x independent with less than 20% fluctuations. Moreover, within 



the same accuracy it is A-independent function as well. Recalling the defenition (|2.11| ) of the 
saturation scale, we obtain 



Qs{x) = Qso{A)x' 



Qn 



0.32 ±0.05. 



(3.14) 



The power is similar to the power q obtained for the proton. Note that the A dependence 
of the saturation scale is found to be x independent: Qso{A) ~ A^'^, where pi is a constant. As 
was explained for the proton case, the initial values Qso{A) and hence the power pi cannot be 
deduced from the scaling behavior only. 

The saturation scales according to the definition ( |2.13| ) are shown in the Fig. (^,a) for four nuclei 
Au, Mo, Ca, and Ne. The numerical errors do not exceed 10%. For the sake of comparison 
we present in the Fig. (P,b) the saturation scale obtained from ( p.lO|) |2l|| . From the saturation 
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Figure 6: The saturation scale Qs is plotted as a function of x for four nuclei Au, Mo, Ca, and Ne. 
(a) - The result is obtained from ( 2.lS^ ); (h) - the result is obtained from (2.1C) [21]. 



scale obtained we can deduce its dependence on the atomic number A, where the power law 
Qs ~ jg assumcd. All the nuclei are divided on two groups: light nuclei (A^e, Ca, Zn) and 

heavy nuclei {Zn, Mo, Nd, and Au). The table (|l]) presents the powers p2 for various values of 
X. The relative errors in the table are estimated not to exceed 20%. On one hand, the power p2 
is seen to decrease with decreasing x. This observation agrees with the results of the Ref. [pT| . 
On the other hand, it can be deduced from the table (|^) that within the errors the power ^2(2;) 
can be viewed as a constant and its average value is in a perfect agreement with the value 2/9 - 
the result of the Ref. 
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The scaling phenomena described above reveals themselves in the energetic gain for performing 
experiments on heavy nuclei. The solution for one nuclei at given x coincides with the solution 
for another nuclei but at different x: 



iVA,(rx,x) ~ NAAr^,\{Ai,A2)x) . 



(3.15) 
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The coefficient A occurs to be x independent (for example X{Ne,Au) ^ 5). The relation ( p.l5 ) 
is a consequence of the scaling phenomena implying Qs,Ai{x) ~ Qs^a2{^{^i^ ^2) x). This leads to 
the relation 




(3.16) 



For Au and Ne this relation gives A(A^e, Au) ~ 5 in total agreement with the direct analysis of 
the solutions. 



4 Conclusions 

In the present paper the scaling phenomena in DIS were studied. The research concentrated on 
the nonlinear evolution equation ( |1 . 1|) governing the dynamics. The solutions to this equation 
were recently found numerically in the Ref. []T^ for the proton target and in the Ref. [^] for the 
nuclei. 

A criteria for the scaling based on the solution of the nonlinear evolution was defined and checked 
numerically. For the proton we found scaling in all kinematic regions of study (10-^ <x< 10-2, 
0.25 GeV^ < Q2 < 2.5 x 10^ GeV^) and with a very good accuracy of order few percents (< 5%). 
The result is in agreement with the discovery of the scaling in the experimental data on the 
structure function F2 [jl[. The scaling behavior is predicted to improve at the LHC and THERA 
energies. It is important to note that the scaling phenomena exist also at distances much shorter 
than the saturation scale At very short distances, where the linear evolution takes place, no 
scaling should be observed. Nevertheless we found that this scaling do exist numerically with 
about 10% accuracy and we are not able to detect its violation. 



The solution found in of the nonlinear equation was used to estimate the saturation scale 
Qs{x). In the present work we gave two new definitions of the saturation scale based on the 
scaling phenomena. In spite of considerable uncertainty in the value of the saturation scale all the 
definitions predict that it grows with decreasing x in accordance with the theoretical expectations 
0) 0) H) m • If we allow ourselves to average between all the result for the saturation scale depicted 
on the figure (^ we would obtain a prediction shown in the Fig. (0,a). The relative errors for the 
latter are roughly 30% for all x which indicate the uncertainty in the saturation scale definition. 

The scaling phenomena were observed for the nuclear targets, which conffim conclusions of the 
Ref. |^0|. The saturation scale estimated from the scaling displays the power law dependence on 



both the atomic number A and the energy variable x: Qs,Aix) ~ A^^ x . The value obtained 
for the power p2 is in a good agreement with the value 2/9 deduced in the Ref. |20[| . 

Both the values of the nucleus saturation scales and their A dependence obtained in the present 
paper are slightly different from the ones found in the Ref. ^|. The main source of this effect is 



certainly comes from the difference in the saturation scale definitions. Since we do not know what 
definition is better we combine all the information we have, proceeding similarly to the proton 
case. The results of this procedure are presented in the Fig. (0,b). We hope that our predictions 
of the saturation scales Qs,a for various nuclei will serve as a theoretical ground for the RHIC 
data analysis in high parton density QCD 
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X 



(a) (b) 

Figure 7: The average saturation scale Qg as a function of x. (a) - the result for the proton; (h) - the 
result for the nuclei. 
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